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O 

f^ ' Starting from the quantum kinetic equation for the non-condensate atoms and the generalized 

("^ , Gross-Pitaevskii equation for the condensate, we derive the two-fluid hydrodynamic equations of 

^S| . a trapped Bose gas at finite temperatures. We follow the standard Chapman-Enskog procedure, 

starting from a solution of the kinetic equation corresponding to the complete local equilibrium 
between the condensate and the non-condensate components. Our hydrodynamic equations are 
r^ • shown to reduce to a form identical to the well-known Landau-Khalatnikov two-fluid equations, 

with hydrodynamic damping due to the deviation from local equilibrium. The deviation from local 
equilibrium within the thermal cloud gives rise to dissipation associated with shear viscosity and 
thermal conduction. In addition, we show that eff'ects due to the deviation from the diffusive lo- 
cal equilibrium between the condensate and the non-condensate (recently considered by Zaremba, 
r^ ■ Nikuni and Griffin) can be described by four frequency-dependent second viscosity transport coef- 

O ' ficients. We also derive explicit formulas for all the transport coefficients. These results are used to 

introduce two new characteristic relaxation times associated with hydrodynamic damping. These 
relaxation times give the rate at which local equilibrium is reached and hence determine whether 
one is in the two-fluid hydrodynamic region. 
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^ ■ I. INTRODUCTION 

5 \ At very low temperatures, the dynamics of a trapped Bose gas is described by the time-dependent Gross-Pitaevskii 

(GP) equation for the macroscopic wavefunction of the condensate. As discussed in several recent reviews [03], there 
is excellent agreement (within a few percent) between experimental observations of collective modes for T ^O.dTBEC 
and theoretical calculations based on the T = GP equation. At elevated temperatures where the condensate 
is appreciably depleted by thermal excitations, one must consider the coupled motion of the condensate and non- 
condensate degrees of freedom. In a recent paper, Zaremba, Nikuni and GrifBn g] derived a generalized Gross- 
ly ' Pitaevskii equation for the condensate atoms and a quantum kinetic equation for the non-condensate atoms, which 
r^ \ can be used to discuss the coupled dynamics of the two components at finite temperatures. These two components 
^\ . are coupled through mean-field interactions as well as collisions between the atoms. 

^^ ■ Two limiting cases for the dynamics of the gas correspond to the collisionless and hydrodynamic regimes |g,^ . Up to 

^^ , the present, most experiments on the collective modes of Bose gases are thought to be in the low-density collisionless 

^^ limit. In this regime, the main effect of the non-condensate (thermal cloud) component is to damp the condensate 

oscillations. In addition to Landau and Beliaev damping due to the dynamic mean-field interaction between two 

components, there is damping arising from the lack of diffusive local equilibrium between the condensate and the 

C non-condensate H . The latter mechanism of damping has recently been worked out in detail by Williams and Griffin 

^ §•.... 

(~| In contrast, in the regime where collisions between atoms arc rapid enough to establish a state of dynamic local 

O • equilibrium in the non-condensate gas, the dynamics of the system is described by hydrodynamic equations for a 

O ' few local variables. This regime contains much new physics and it should be accessible taking advantage of the 

J> \ larger densities now available as well as the large scattering cross sections close to a Feshbach resonance B. We 

r '~j I have recently given a detailed derivation and discussion of two-fluid hydrodynamic equations for trapped atomic 

rN \ gases at finite temperatures P,p|,Kp|,p^ , starting from a generalized GP equation coupled with a kinetic equation. 

S I These equations were derived under the assumption that the non-condensate atoms are in local thermodynamic 

" " ' equilibrium among themselves but are not in diffusive equilibrium with the condensate atoms. The resulting ZGN' 

hydrodynamic equations P,p[ involve a new characteristic relaxation time t^, which is the time scale on which local 

diffusive equilibrium is established. This equilibration process leads to a novel damping mechanism which is associated 

with the coUisional exchange of atoms between the two components. This ZGN' hydrodynamics is briefly reviewed in 

Section ||. In Ref . [|lO[ , we generalized the ZGN' equations to include the effect of deviations from local equilibrium, 

and worked out hydrodynamic damping associated with the collisions among the non-condensate atoms. At finite 

temperatures of interest, this deviation from local equilibrium within the thermal cloud gives rise to damping associated 

with thermal conductivity and the shear viscosity. Such a generalization was first discussed in Section V of Ref. IIj] 



starting from the ZGN hydrodynamic equations. We also note that the thermal conductivity and shear viscosity 
were first derived for a uniform Bose-condensed gas at finite temperatures in pioneering papers by Kirkpatrick and 
Dorfman [l^ ]. 

In the present paper, building on our recent work with Zaremba, we present a more complete derivation of two- 



fluid hydrodynamic equations, including dissipation. In Section III, we solve the kinetic equation by expanding 
the distribution function /(r,p,t) around /^"^(r, p, i), the latter describing complete local equilibrium between the 
condensate and the non-condensate. We follow the standard Chapman-Enskog method used to derive hydrodynamic 
damping in the kinetic theory of classical gases. In this treatment, the lowest order hydrodynamic equations involve 
no dissipative terms. All hydrodynamic damping effects are included by taking into account deviations from the local 



equilibrium distribution f^*^'. In Section IV, we prove that, with appropriate definitions of various thermodynamic 
variables, our two- fluid hydrodynamic equations with damping have precisely the structure of those first derived 
by Landau and Khalatnikov ||ll|,|T2[. In particular, the damping associated with the collisional exchange of atoms 
between the two components, which has been discussed at length in our previous work [pp|JlO||, is now expressed in 
terms of frequency- dependent second viscosity coefficients. This type of damping is in addition to the usual kind of 
hydrodynamic damping associated with shear viscosity and thermal conductivity |lC| , |l4| . 

In Section M, we also derive explicit expressions for all the transport coefficients that appear in the dissipative terms 
in our two-fluid hydrodynamic equations. For the purpose of illustration, we evaluate the temperature dependence 
of all the transport coefficients for the case of a uniform Bose gas. The case of a trapped Bose gas is quite different. 
Because the condensate density is always much larger than the non-condensate density in the central regions of a 
trapped Bose gas, we find that the collisions between the condensate and non-condensate atoms are the dominant 
contribution to all transport coefficients. 

The lengthy analysis in Sections [II, IV and M is, of necessity, very complex and conceptionally quite subtle. 
Most readers will only be interested in the final conclusions, which we now summarize. We prove that the two-fluid 
hydrodynamic equations of a trapped Bose gas can be written precisely in the well-known Landau-Khalatnikov form, 
as summarized in (p6|). The new feature which arises in a trapped Bose gas (as compared to superfluid ^He) is 



that, as noted above, the four second viscosity coefficients can b e fre q uenc y-depe nden t, as given by ( 106 ). Finally, 



explicit expressions for all the transport coefficients are given in (II£), (f3C) and (139). An important final result of 
our analysis (see Section 0) is a precise definition of three characteristic relaxation times associated with the various 
transport coefficients. Moreover we show that the collisions between the condensate and non-condensate atoms always 
play the dominant role in the hydrodynamic damping of trapped Bose gases. 

The present paper bririgs to a natural conclusion our series of papers (with Zaremba) on the two-fluid hydrodynamics 
of trapped Bose gases pj9| |l(]| , [l4| . Much remains to be done solving these hydrodynamic equations and experimentally 
checking the predictions. 

II. ZGN' HYDRODYNAMIC EQUATIONS 

In this section, we derive the most general form of hydrodynamic equations for the condensate and non-condensate 
at flnite temperatures, and then briefly review the ZGN' two-fluid hydrodynamics. We start with the underlying 
coupled dynamical equations for the non-condensate and the condensate, as derived in Ref. and recently reviewed 
in Ref. Q. The non-condensate atoms are described by the distribution function /(r, p, t), which obeys the quantum 
kinetic equation 

9/(r,p,t) p ,. ^ .. ^ ,, ,. 

T77 + — • Vr/ r, p, t) - Vrt/ ' Vp/ r, p, t) 

mm 

= Ci2[/,<i>]+C22[/]. (1) 

Here the effective potential U{r,t) = C/ext(r) + 2g[nc{j^,t) + n{r,t)] includes the self-consistent Hartree-Fock (HF) 
mean field, and as usual, we treat the inter-atomic interaction in the s-wave approximation with g = Airh a/m. The 
condensate density is 7ic(r,t) = |$(r, i)p and the non-condensate density n{r,t) is given by 

^(^'^)=/(W^(^'^'^)- (2) 

The two collision terms in (|l|) are given by 

CHf]^-^ J dp, J dp, J dp. 
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(4) 

with / = /(r, p, t), fi = /(r, Pj, i). The expression in ^ takes into account the fact that a condensate atom locally 

has energy ec(r, i) = /ic(r, i) + i2'mv'^{r,t) and momentum rnVc, where the condensate chemical potential Hc and 

velocity Vc will be defined shortly. On the other hand, in our finite-temperature model, a non-condensate atom 

2 
locally has the HF energy ep{r,t) — 2:^ + U{r,t). This particle- like dispersion relation limits our entire analysis to 

finite temperatures. 

The equation of motion for the condensate is given by a generalized Gross-Pitaevskii equation [pj for the macroscopic 

wavefunction <E>(r,t) (see also Ref. [[15[) 
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C/cxt(r) + gndr, t) + 2gfi(v, t) - iR{r, t) 



1>(r,t), 
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rCi2[/(r,p,t),<J>(r,i)]. 



(5) 



(6) 



(7) 



The dissipative term R in (5) is associated with the exchange of atoms between the condensate and non-condensate, 
as described by the collision integral Ci2[/, $] in (4). We see that (1) and (5) must be solved self-consistently. It is 
customary to rewrite the GP equation (5) in terms of the amplitude and phase of $(r, i) — ■y/nc(r, f)e'^*-'"'*', which 
leads to (vc = KV9{r,t)/m) 



~dt 
d 



V ■ (n^Vc 



-ri2 [/,<&], 



"^ ( ^ + Vc • V ) Vc = - V^c 



where the condensate chemical potential is given by 



Aic(r,t) 



2my^nc{r,t) 



+ C/oxt(r) + gnc{r, t) + 2gn{v, t) 



(8a) 
(8b) 

(9) 



One sees that ri2 in (8) plays the role of a "source function" in the continuity equation for the condensate, arising 
from the fact that C12 collisions do not conserve the number of condensate atoms ||] . 

We note that the set of equations (l])-(R) has also been derived using the elegant Kadanoff-Baym formulation 
|L6|-[18| . More recently, this KB derivation has been extended to cover low temperatures as well [|l^ , by working with 
a Bogoliubov quasiparticle spectrum instead of the simpler HF spectrum used in the present paper. One could use 
this extension as the basis for generalizing the present paper. 

Following the standard procedure in the classical kinetic theory of gases pQ] , we take moments of the kinetic equation 
(|l|) with respect to l,p and p^ to derive the most general form of "hydrodynamic equations" for the non-condensate. 
These moment equations take the form (/i and v are Cartesian components): 



— + V-(nv„) = ri2[/], 

mn I — + v„ • V I w„^ = — n- m(w„^ - Vc^Ji)Tl2[J\ , 

'1 



de 
di 



V Q D.,P., 



l(v„ - Vc) + ^ic-U 



ri2[/]. 



(10a) 
(10b) 

(10c) 



The non-condensate density was defined earlier in (yj), while the non-condensate local velocity is defined by 

-(MK(M)./^^/(r,p,t). (11) 

In addition, we have introduced the following quantities: 



Q(r,t)^ / _^ (p_^v„)2 ii_v„ /(r,p,t), (12b) 



dp 1 , n2 / P 

(ZTTn)-^ 2m \m 



Finally, the symmetric rate-of-strain tensor appearing in (10) is defined as 



The hydrodynamic equations in ( |lO| ) are formally exact, but obviously are not closed. To proceed, we must specify 
the conditions under which the dynamics of the system to be determined. In the ZGN' hydrodynamics, we consider 
the situation in which the C22 collisions are sufficiently rapid to establish local equilibrium among the non-condensate 
atoms. This situation is described by the local equilibrium Bose distribution for the thermal cloud, 

/>^P^O= ^,t^(p_^v„W,]_^ - (14) 

Here, the temperature parameter /3, normal fluid velocity v„, chemical potential fl, and mean field U are all functions 
of r and t. One may immediately verify that / has precisely the form such that it satisfies the condition C22[f] = 0, 
independent of the value of jl. In contrast, one finds that Ci2[f] does not vanish in general, namely 

ajf] = J^!!!^ri _e-/3(A-i™(v„-v.)^-Mc)i 

''^■'^ (27r)2;i4[^ J 

X dpi dp2 / dp3(5(mvc + Pi - P2 - P3)S{^i + ^c - £2 - £3) 

X [^(p - Pi) - S{p - P2) - <5(p - P3)](l + /i)/2/3. (15) 

Using the local distribution function (14) to evaluate the moments in (2) and (12), we find that the heat current 
Q(r,t) =0, and that 

f dp ~ 1 

n{r, t)^ J .^^^.3 /(r, p, t) = ^53/2(2) , (16) 

P,.{r,t)=S,.P{r,t) ^ 5,,j ^_|^lP_^/(r,p,i) = J^. — 95/2(z). (17) 

Here z{v,t) = e'3['i('-'*)-'^('-^*)l is the local fugacity, A(r,t) = [2T:h'^ /mkBT{r,t)]^/'^ is the local thermal de Broglie 
wavelength and gn{z) — J^'iZi z''/l" are the Bose-Einstcin functions. The kinetic energy density is given by e(r, i) = 
|P(r,t). 

To summarize, using / ~ /, we obtain the ZGN' lowest-order hydrodynamic equations for the non-condensate given 
in Refs. @,H,P[. In the linearized version of these ZGN' hydrodynamic equations, the condensate is described by 

/Inn 

i^ -f V • (n,o<5ve) = -Sr,2 , (18a) 

m-^ = -V5fi,, (18b) 



and the non-condensate variables satisfy 

or — 

— J^ + V • (no5v„) = (5ri2 , (19a) 

at 

or 

mhn—^ = -VSP - SfiVUo - 2gnoV{5n + Sric) , (19b) 

dSP 5 ~ 2 ~ 2 

^ _ V • {PoSvn) + -v„ • VPo - ^ffncO'5ri2. (19c) 

ai 3 3 3 

The fluctuation of the condensate chemical potential is given by 

(5/ic = gSuc + 2gSh. (20) 

This assumes the Thomas- Fermi approximation, which means that the first term in (m, the quantum pressure term, 
is neglected. The source function i5ri2 can be usefully expressed Q in terms of the difference between the condensate 
and non-condensate chemical potentials, namely 

6Ti2^ —{Sfj.~Snc) = ^(^Mdiff, (21) 

Tl2 Tl2 

where fi^is ^ P- ~ fJ-c and ri2 is the mean collision time MM associated with C12: 

T12 (27r)5r J J J 

X(5(pi - P2 - P3)S{flcO + £1 - £2 - £3) 

x(l + /io)/2o/30, (22) 

where fio = \^e^o{£i-i-iaa) _ j^j-i jg ^j^g static equilibrium distribution function and ii — ^ -I- Uq. 

Since one can show P,^Q] that (5/idiff can be written as a linear combination of 6n and 5P, the above coupled 
hydrodynamic equations in (|8|) and ( [l9| ) are a closed set for the variables 6nc,Svc,5h,Svn and SP. However, it will 
be useful later to have an equation of motion for (5/idiff- This is given by [see Eq. (86) of Ref. Isj] 

— tt; — = -g V • nco(OVc - dv„) - -grico V • ()v„ . (23) 

ot 3 T^ 

Here r^ is a new relaxation time governing how (5/idiff relaxes to zero, i.e., how fast diffusive equilibrium is reached 
between the condensate and non-condensate components. It is related to the collision time T12 in (E3) by the expression 

1 Pgr^^ (24) 

r^(r) (THT12 

where the dimensionless hydrodynamic renormalization factor (Th is given by 



f^Hir) 



|Po -I- 2gnonco + |g^7oneO 



1 



-1 



(25) 



i^o7o - fg^o 

where 7= j^§^9i/2iz)- 

We can now look for normal mode solutions of the linearized ZGN' equations in ([iq ) and (p!9|). Assuming that these 
fluctuations have a time dependence e~*'^*, one can solve ( p3| ) for (5/idiff to give 

<5Aidiff = --; \ .gV • [nco((5vc - (5v„)] + -gucoV ■ (5v„ } . (26) 

1 - lUJT^ [ S ) 

In the limit lut^ —f 0, one sees that Sfidie —^ 0. This situation corresponds to the complete local equilibrium between 
the condensate and non-condensate components, with /i(r,t) = /Xc(r, i). In this limit, one can prove that our ZGN' 
hydrodynamics reduces to the Landau two-fluid hydrodynamics without dissipation terms, as discussed in detail in 
Ref. 1^ . It is clear that fluctuation of Sfidis described by ( p3| ) [or equivalently (^)] gives rise to a new relaxational 
mode in addition to usual collective oscillations of the condensate and non-condensate (for a uniform superfluid. 



these are the first and second sound modes). For a uniform Bose gas, the frequency of this new mode is given by 
wr = —iJTy, [|9|]. Thus, in general, our ZGN' equations predict the existence of a new relaxational mode, associated 
with the equilibration of the condensate and non-condensate collective variables. 

In Ref. pel , we have extended the theory to include small deviations from the local equilibrium distribution / in 
(14) . These give rise to new dissipative terms in the non-condensate equations associated with the shear viscosity (77) 
and the thermal conductivity (k) of the thermal cloud. The damping of first sound, second sound and the relaxational 
mode due to the effect of normal fluid transport coefficients was calculated in Ref. [Q. In particular, we showed there 
that the relaxational mode was strongly coupled to (and renormalized by) fluctuations in the local temperature and 
hence the thermal conductivity. 

In summary, the ZGN' hydrodynamics exhibit the physics of the coupled dynamics of the condensate and non- 
condensate atoms in a clear fashion. However the approach used in the ZGN' theory has a disadvantage in that it 
is not based on a small expansion parameter, in contrast to the more systematic Chapman-Enskog procedure used 
here. In Ref. Eoj, we only included the effect of C22 collisions to discussing the deviation from local equilibrium. 
This neglect of C12 collisions in this connection is only justified when the condensate density is very small compared 
to the non-condensate density (since the C12 term in ^^ is proportional to ric). However in a trapped gas, the C12 
collision integral is always significant since the condensate is strongly peaked at the trap center, with a density much 
larger than the non-condensate even at temperatures close to Tbec- Thus, we must treat both C12 and C22 when 
considering deviations from local equilibrium. 

In the following section, we present a more systematic derivation of the two-fluid hydrodynamics, by following the 
standard Chapman-Enskog procedure. This derivation is similar to the work by Kirkpatrick and Dorfman |l3[ for a 



uniform Bose gas. As we discuss in Section [V, this new approach allows us to show that the extended ZGNM;heory 
can be written in a form completely equivalent to the Landau-Khalatnikov two-fluid hydrodynamics |l2] when we 
include hydrodynamic damping. This set of equations involves the thermal conductivity, shear viscosity and four 
frequency-dependent second viscosity coefficients. The latter are shown to arise from the fact that the condensate 
and non-condensate are not in diffusive equilibrium (^c 7^ A)j as discussed by ZGN' j|Jg]. 

III. CHAPMAN-ENSKOG EXPANSION FOR A BOSE-CONDENSED GAS 

A. Low^est-order hydrodynamic equations 

Following the standard procedure of the Chapman-Enskog expansion |24| ] , we introduce a small expansion parameter 
a and rewrite the kinetic equation dl^) as 



df{r,p,t) ^ p 
m 



dt +f •Vr/(r,p,i)-VrC/-Vp/(r,p,t) 



-(Cl2[/]+C22[/]). (27) 

a 



This expansion parameter a will be eventually taken to be 1, but allows one to develop a perturbative solution of 
(pTI). In order to solve the quantum kinetic equation, we formally expand the distribution function /(r, p, t) in powers 
of a: 

f = /^"^ + af^'^ + ■■■. (28) 

Using this expansion (Eq), we can also expand the various hydrodynamic variables in (nfl) 

n = n(0)+«nW + ..., P^, = pW + «pW + . . . ^ q ^ q(o) + ^q(i) + . . . ^ 

i = f:(°) + a?(i) + • • • . (29) 

The superscript (0) denotes the local equilibrium solution (see below) which is determined by the collision integrals 
(formally when a ^ 0). We also redefine the source function ri2 in (0) as 

- /■ ^^Ci2[/, $] ^ i(r(°^ + «rW + a'T,, ■ ■ .). (30) 



We also have an expansion for the condensate wavefunction 

$ = $(")+a$(i) + ---. (31) 

6 



In this expansion, however, we assume that the total local density n{= ric + n) is not altered by the higher order 
correction terms /^*' {i > 1) in (|2q). That is, we have 



Tic = nf^ + an(i) + • • • , (32) 

but 

n = 4") + n(°), n« + n(*) = (i > 1). (33) 

We also assume that non-local correction terms /'•*■' make no contribution to the non-condensate velocity fields v„ 
or to the phase 9 of the condensate wavcfunction (and hence to the condensate velocity Vc). Finally, the condensate 
chemical potential in (S) (we work within the Thomas-Fermi approximation) is given in the expansion 

^lc{r, t) = L/c,t(r) + .g[n(r, t) + n(r, t)] = ^l';^\r, t) + a^(i) (r, t) -f • • • , (34) 

with 

A.(") = f/ext + 9{n + n(")), ^i^^ = gh^'\ (35) 

Using the expansion ( p8| ) in the kinetic equation (|2^), we find that the lowest order solution /(°) is determined from 

Ci2[/(o),$(o)]-fC22[/(°)]=0. (36) 

The unique solution of (^^ is given by the "diffusive local equilibrium" Bose distribution function, namely 

7 (r, P, i) - ^^(r,t)[_i_(p_™v„(r,t))2 + C/(r,t)-/i(0)(r,t)] _ ^ ' ^^'> 

Here the local equilibrium non-condensate chemical potential /i'^*'^ is given by the condition that Ci2[/*-°-', <&*-°^] — 0, 
which gives 

A("^=/^i"^ + f(v„-v,)^ (38) 

Using (psj), this is equivalent to 

^(0) = ;.(") -,- ^(v„ - v,)2 = C/e.t + 5n + 5n(°) + y (v„ - v,)^ (39) 



in conjunction with 



^'''(-'^)- I J^f'"'i-^P^*)-T^9V^(^'"'^- (40) 



Here z^°\r,t) = e^^^ (i",t)-C'(r,t)] jg ^^j^g local fugacity in diffusive local equilibrium. 

It is important to appreciate that diffusive equilibrium is not defined by the distribution function f^'~'' alone, but 
is determined self-consistently with the non-condensate chemical potential as given by (|38|) . One may immediately 
verify that /(°) satisfies C22 [/'"'] = 0, independent of the value of p,^°\ In contrast, C\2[f^°\^^°^ = only if 
the local chemical potential of the thermal cloud is given by ( pq ) and the condensate and non-condensate densities 

are determined self-consistently. Of course, it immediately follows that since Ci2[/^*'-', <I'^°-'] = 0, we have TI2 — 
Ti2[f'-°\^'-°^] = and hence (|o|) reduces to 

ri2 = riy + arf^^ + • • • . (4i) 



Using the local distribution function (pTJ ) to evaluate the moments in ( |l2b| ), we find that the heat current Q^"^ (r, t) = 
0, and 

P(°;(r,t) =<5,.p(")(r,t) ^ V/ (^^^^^/'"Hr,p,t) = <5,.^g5/2(.('')). (42) 

The local kinetic energy density is given by e(°^(r, t) ~ |p("^(r, t). 



To summarize, the lowest-order hydrodynamic equations for the non-condensate are given by 



f + V.(nv„)=rW, 



d_ 
di 



V v„ 



dP 
'dt 



V-(Pv„) = -^PV-v„ + ^ 



VP - nVU - m(v„ - Vc)r 
1 



(1) 

12 I 



-m{Wn - Vc) + ^Jic"U 



.(1) 



(43a) 
(43b) 

(43c) 



where n = n(°\ P = p(°) and fie — fJ-c are given by (|3q), ( [40| ) and (42). It should be noted that the above equations 



involve the source term Fj^j • Even though Ci2[/'-°''] = 0, one sees from (41) that the lowest order contribution is in 



-(1) 



fact given by ^^2^ which involves the contribution from the next order correction f^^'. Later we will derive an explicit 
expression for ri2 when we include the effect of deviations from the local equilibrium distribution and transport 
processes. Here we only display the result for the lowest order contribution which enter into ([43[) [see also 



r('2^(r,i)=f7ff <j V-K(Ve-V„)] 



1 



(44) 



.(1) 



the 



where an is defined by (BE 

It is important to note that even though r^^ involves an integral over the collision integral C12 (see (|3C 
expression in ( p4| ) does not involve any collision time. The expression for T]^2 iii (44) is consistent with the ZGN' 
result for SfidiS given in (|26| ) in the limit wr^ -^ [using ( |2l| ) and (|2^)]. We recall that in this limit, one has S^diS -^ 
and thus / in (Ol) reduces to /'■°-' in (p^). Therefore the hydrodynamic equations given in (E^) are equivalent, in the 
i^Tfj^ — > limit, to those given by the ZGN' theory. As noted in Ref. [^, these coupled lowest-order hydrodynamic 
equations in (Q) and (H) can be combined and also shown to be precisely equivalent to the Landau two-fluid equations 
without dissipation due to the transport processes. In Section |rv| , we prove this equivalence in the more general case 
when dissipation is included. 



B. Two-fluid equations with hydrodynamic dissipation 



We next consider the deviation (g8|) from the local equilibrium distribution function /'^"^ to first order in the 
Chapman-Enskog expansion. This deviation f^*^' gives rise to additional dissipative terms in the hydrodynamic 
equations. As usual, in determining the dissipative terms, we restrict ourselves to terms of first order in the velocity 
fields v„ and Vc. Following Refs. [|l^,|lj,|l^ , we write the first correction term in ( |2q ) in the form 



/W = /(")(r, p, i)[l + /W(r, p, i)]V'(r, p, i), 
and work with 'ip{r, p, t). To first order in a, the C22 and C12 collision terms in 



(45) 



reduce to {f ^ f(°^ + af^^'^). 



-C22[/]^C22[/(^'] 



29' 



{2TTfn 



5%r 



dp2 / dps / dp4 



-£p4) 



X(5(p -K P2 - P3 - P4)'5(epi + ip^ ~ ep3 

X/(«)/f (1 + /f ) (1 + /f ) (^3 + V4 - ^2 - ^) 



(46) 



-Cx2\!M 

a 



(27r)2?i4 
X(5(mVc 



dpi / dp2 / dp. 



Pi -P2 -P3)Sie 



(0) 



X [S{p -pi)-S{p- P2) - S{p - Pa)] 

x(l + /}"Vi°Vf (-/3A^i'^ + V'2 + ^^3 - ^^l) 
= -/55n(i)Li2[l]+Li2[^], 



(47) 



where Ec — iJ-c 



M 



rav"^ and ^c is given by (35). The linearized L12 operator is defined by 



^^^[^]"(i^/'p^/'p^/'p^ 



xS{m-Vc + Pi - P2 - P3)^(ei°^ + £pi - £p2 - Sps) 
X [5{p - pi) - S{p - P2) - ^(p - Pa)] 

x(l + /l"^)/f/f (V'2+V'3-V'i). (48) 



Using ([46|)-(|48|) and expanding the kinetic equation (|27|) to first order in a, we find that the first non-local correction 
f^^' is determined by the equation 

^"^''P'^^ + P . Vr/(°) (r, p, t) - VrC/ • Vp/(°) (r, p, t) 
at m 

= -l3gn^^^L^2[l] + ^12^ + ^22[V']- (49) 



Here d'~' /dt means that we use the lowest order hydrodynamic equations given by MS) in evaluating time derivatives 
of v„, /i, T and U. The resulting linearized equation which determines the function tp is (for details, see Appendix A) 



u- VT 

T 



mu'^ 535/2(2) 



2kBT 233/2(2) 



kBT 



Dfj^u UpU^ - -(5, 



/it/ <- 



;ai) ^ I /("Hi + /^°^) + /35"^'^ii2[l] = il2M + i22M ^ iM, 



^2 + ^^ai ) ^ } /^^'Hl + /^'O + /35"^'^ii2[l] = il2M + i22M = iM, (50) 

where the thermal velocity u is defined by mu = p — ?tiv„ and z = z^^' . The dimensionless thermodynamic functions 
(Ti, (72 in JSOJ ) are defined by 

^^(0)^(0)[^(0)_[;]„|[^(0)]2 

'^^'■'^'*) = |p(0)^(0)_|[^(0)]2 ' 

5pfo)^(0) _ [^(o')]2 r~(0) „ [;! 

--(-'^^-^ iP(0),(l)-|M°)]- ' ^''^ 

where 7^°)(r,i) = ^3i/2(z'-"-'(r, t)) = T*-"-*/.?- We note that C12 enters in three separate places in (pQ). 
The linearized collision operators L12 and ^22 satisfy the conditions 

Li2[p-mv,] = 0, Li2[£p-ei")]-0, 

i22[l]=0, L22[p]=0, L22[ep]^0. (52) 



In order to have a unique solution of ( |5C| ) for ■0, we impose the following additional constraints: 

y'dpu/(")(l + /("))^ = 0, (53a) 

y'dp(^u2 + C/-/l("))/(")(l + /("))^ 

= iy'dpln(l + /(")"')/(!) =0. (53b) 

Physically, the constraint (|53a| ) means that the deviation from local equilibrium make no con tribution to the local 
velocity field v„ defined in (|1 1|) . As we discuss in more detail in Section |^, the constraint (53b) means that the total 



energy density and the local entropy density are not altered by the deviation f^^'. They have the the same value as 
given by r^' . 

Since ( pO| ) is a linear integral equation for ip, one may write the most general solution in the following form EJ]: 

^^•"- ^ - ^ 1,2, U,„^^U1), 



V'(r, p, t) = -^^A{u) + D^, \u^u, - -^u'S^u] B{u) + T'^^D{u), (54) 

where the dependence on (r, t) is left implicit and u^ is a component of the thermal velocity. Here the functions 
A{u), B{u) and D{u) are given by the solution to the following three linearized integral equations: 



TOW^ 5^5/2(2) 



/(")(! + /(")) = L[uA(u)], 



2kBT 253/2 (z) 



TO 



For the constraints (m) to be satisfied, we also need to require that 

dp 



{2Trhf'' '^ ' ■' ' \ 2 ' J 



{2-Kny 



rK^ + fn 



(0)^ f ™^^^^_ ~(0)\ 



+ U - fl'-"> D{u) = 0. 



(55a) 
(55b) 
(55c) 

(56a) 
(56b) 



Using the solution for ip given in (p4[), one finds that the corrections due to /'^' in ( [45[ ) to the various hydrodynamic 
variables are given by 



n(i) 



'^P /(O)(i + /(o))^(,)r«(r,0^-rr«, 



{2TTh) 






£'a"' ^ -TrDfJ^^ 



with 



P«=r^(L/-A(°))rW^r^5n?M2\ 

g(i) ^ 3^(1) 
2 



-(1) 



(57) 

(58) 
(59) 

(60) 
(61) 



We note that n and P are both altered by an amount proportional to Fj^j • The transport coefficients ry and k are 
associated with the functions A{u) and B{u), 



V 



dp 



f(o) 



(0)^ 



157 (27r;i)^ 
dp 



6T7 (27r?i)3 

The relaxation time r defined in (p7f), namely 

dp 



«4A(u)/(")(i +/(«)). 



(27r;i)' 



-/(«)(l + /("))i^(u), 



(62a) 
(62b) 

(63) 



plays a crucial role in the subsequent analysis. Using ( pTj ) in (p5q), one can rewrite the integral equation for D{u) in 
the form 



0-2 



UbT V n(o) 



1 



Ko) 



'Ti ^7m /'"'(I + /^'") - r/35^i2[l] - L[i?(«) 



(64) 



In Section M, we solve the three linearized equations listed in ( [55| ) . It will be shown there that the solution for the 
function DTu) is 



D{u)^-^y<r2 , 3^^^ 



O"! 



(65) 



Using this, one finds that r can be identified with the relaxation time r^ defined in (P4|). In the present discussion, 

the physical meaning of the relaxation time r^ can be clearly seen by writing the source function F-j^j iii the form (see 
(§) and (H)) 



10 



il2 --^ 



n(0) 



(66) 



This kind of relaxation term in the two-fluid hydrodynaniic equations such as ( pa| ) and (43a) was also discussed in 
a pioneering paper by Miyake and Yamada |22| ] in discussing the liquid "^He near the superfluid transition (where a 
phenomenological relaxation time equivalent to r^ was introduced) . 

In summary, we have obtained the following hydrodynamic equations for the non-condensate including the normal 
fluid transport coefficients (we now set the expansion parameter a = 1): 



dn 



V • (nv„) = ri2 



■mn\ ^ + v„ • V j w„^ 
d 



dP dU 

1^ + '^1^ = -m(w„^ 



2r] 



de 
dt 



+ V • (ev„) 



dxi, 

(V • v„)P = 






Vc^i)T 



'cfj.)^ 12 



"i(v„ - Vc)^ + ^c-U 



+V ■ (kVT) + 2t] 



ri2 

2 



D^,--{TtD)5^, 



where h and P are given by 



P = P^'^+r,lgn(M\\ 



(67a) 



(67b) 



(67c) 



(68) 
(69) 



and e = |P. Here n^^^ and P*-"' are given by (|4C|) and (|42|), respectively. The equivalent "quantum" hydrodynamic 
equations for the condensate are given in (H) , where the condensate chemical potential is given by 



,(0) 






We now derive an equation for the function ri2. Using (|68|) and ( |67a| ), one obtains 

5ri2 






dt " dt 

Using the explicit expression for fv-^i given in (HO), one obtains 



V ■{nWn)+Ti2. 



1 dT ( 

2- + 57nc|-^ + o^( 



dfi--^^ dn 



91 



dt \2'^ ' '"■'■7 T dt ' "' \ dt dt 
Using the continuity equation for the total density n, one finds that (^) reduces to 

9^(0) 1 



dt 



Using (|73D in ([TlJ), one finds 

ari2 



dt 



ri2 = 



1-57 



1 



3~ \ 1 c»P „ ,. 

-n + g-fiic 1 y "a^ + " ■ i"'^" + "cVc) 



1 -.97 



3~ \ 1 9r „ ,. 



+ V • (nv„) 



We next use (p9h in the equation for P (given by (67c)) 
dP 5P(o) 2 



dt 



dTi2 



^PV • v„ - ^371 J12 + ^ V • {k,VT). 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



Using the expression for P*^"^ given in (p2|), one finds 
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-^- = l-P + gnnA -- 



gn- 






Substituting this into (^5|) in conjunction with ([7lD, we obtain 



-guc + gn 






gnV ■ (nv„ +ncVc)- 



1 dT 



(76) 



5 ~ 

-P + gfific ] ^-^+ anV ■ (ncVc 

|pV-v„-^V-(acVT). 



VP-v„ 



(77) 



One may now combine (|7J) and (^TJ) to eliminate dT/dt from these two equations. After a certain amount of 
rearrangement, we finaUy obtain our desired equation of motion for ri2: 



dV. 



dt 



ri2 = CTfl i V • [nc(vc - v„ 



1 



-UcV ■ V,; 



2 gggj 
3^ 



V • (kVT). 



(78) 



If we keep the expansion parameter a and expand ri2 as in (|41|), namely ri2 = T\ 



(1) 



-(1) 



we find r-j^2 is given by (}Am and 



ri- 



(2) _ 



-^^m^ 



iy-|^^v.(.VT). 

3 .gn 



arfa^ (we recall that T^"^ = 0), 



(79) 



In closing this section, we discuss the relation between the analysis given in this Section and the ZGN' thcor > 
reviewed in Section ||. In this Section, we started with the complete local equilibrium distribution given by ([3'^. 
We then included the deviation from local equilibrium, as given by (Esh with (p4). We showed that the deviation 
from f^^' associated with D{u) in ( p4[) gives rise to the corrections to the local thermodynamic quantities n, P and e. 
Such corrections did not arise when we included the deviation from / in the ZGN' hydrodynamics [ [L0[ . However, one 
can show that the type of contribution associated with D{u) is, in fact, already contained in the lowest-order ZGN' 
distribution function / given by (|l4|). To see this, it is convenient to linearize the distribution function around static 
equilibrium, using f ~ fo + 5f . In the ZGN' theory [0, one can show that 



(5/ = /?o/o(l + /o) 



5T ( p 



Tq V 2to 



+ t/o - Aico + P • v„ - 2gSn + Sp, 



(80) 



Here we have denoted the temperature fluctuation as ST to make a distinction from the temperature defined in the 
diffusive local equilibrium distribution function (37) (we will find that 6T ^ ST). In the present theory, in contrast, 
one finds (ignoring the terms in (|54| ) associated with the functions A and B) 



Sf = Pofoil + fo) 



ST f p 



To \ 2m 
-fo{l + fo)D{u)STi2- 



+ Uo-^lco] +P-^n- 2gSn + S^f^ 



(81) 



The first term in ( pl| ) represents the deviation from /o included in f^^> while the second term is due to f^^'. Using 
the explicit solution for D{u) given by ( |6q ) (derived in Section 0), we find that ( |8l| ) can be written as 



5f = /3o/o(l + /o) 



ST 2aiTn 



To 3ho 



'5^12 ) ( ^ + ^0 - AicO 



p • v„ — 2gSn 



(82) 



We note that this linearized distribution function has the same form as the ZGN' distribution function in (pO|), but 
with a renormalized local temperature 



Sf = ST~lTo^Sri2, 
3 no 



(83) 



and a renormalized local chemical potential 



12 



Sfl = Sn^°^-giaH' + l)T^STi2. (84) 

Using Sfic — Sfic ' — gT^6Ti2 (see (|70|)) and Sp, from ([8^), we obtain 

S^ditf = 5ji- S^c ^5ji- [(5^^°^ - gTp,5V 12] == ga^^ST 12. (85) 

This relation between <5/idiff and (5ri2 is precisely that given by (^l]) and (p4|), as derived in the ZGN' theory. The 
physical significance of the renormalized thermodynamic quantities, as given by ( p3|) and (p4|), will become clear in 



Section IV 



IV. EQUIVALENCE TO LANDAU-KHALATNIKOV TWO-FLUID EQUATIONS WITH DISSIPATION 



In this Section, we prove that our hydrodynamic equations in (ph and (p7h can be written in the form of the 
Landau-Khalatnikov two-fluid equations. We first display the complete Landau-Khalatnikov two-fluid equations 
involving dissipative terms |l4| : 

^ + V-j = 0, (86a) 

^U d dUext 

TO—— -I- - — [o^iji^ + mnVnf^Vnu + mncVcfj.Vci,) + n 



dt ' dx,.' """ "^ '"^ —'^-'^^ ' Qj.^ 



dXy ' 



D^,,-U^,{TvD) 



<5m.(Ci V • [TOne(ve - v„)] + C2V • v„) L (86b) 



^ = -V (^ + ^ - C3V • [TOn,(v, - v„)] - C4V • v4 , (86c) 

at yni Z ' 

^ + V.(^,sv„-^j.-. (86d) 

The total current is given by j = ric'Vc + nVn and the dissipative function describing the entropy production rate is 
given by [|| 

Rs = C2(V • v„)2 + 2Ci V • v„V • [mricivc ~ v„)] + Cs (V • [mn^iv^ ~ v„)])^ 



+2r] 



1 



-Dp. - 3'5^.(Tri?) 



-(Vr)2. (87) 



As we have discussed in Ref [pi, the normal fluid and the superfluid densities that appear in the standard Landau 
two-fluid theory can be identified with the corresponding non-condensate and condensate densities, within the context 
of our finite temperature model based on the HF approximation for single-particle excitations. We have explicitly 
made use of this correspondence in writing (|8q). We also note that in (|8q ) and (|87|), one can write ncivc — v„) in the 
equivalent form (j — nv„), which is often used. 

The thermodynamic functions that appear in these Landau-Khalatnikov two-fluid equations satisfy the following 
superfluid local thermodynamic relations: 

P + e = fm + sT + mh{vn-Vcf, (88a) 

dP = ndfi + sdT — TO,n(v„ — Vc) • (i(v„ — Vc), (88b) 

de = iidn + Tds + (v„ — Vc) • d[mh(yn — Vc)]. (88c) 

The various local thermodynamic functions which appear in the LK theory have to be carefully defined so that they 
satisfy the relations in (|8^). The local entropy is defined by (as in Ref. 0) 

^[(l + /)ln(l + /)-/ln/]. (89) 



{2-Kn) 
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Using / = /*^°^ + Z*^^) and working to first order in f^^\ one finds 



dp 



{2nh)3 



(1 + /(")) ln(l + /(")) - /(") ln/(") + ln(l + f<-°r')f 



(1) 



(90) 



From the constraint on /'^' given by (53b), one sees that the last term, which arises from /'^', makes no contribution 
to the local entropy. One thus obtains 



1 



5p(0)_~(0)(^(0)_f^) 



1 

T 



2 



n(°) (;.(") 



mn 



(91) 



where we have used n'"' = h + 0{vn, Vc). 

The local energy density e in the Landau-Khalatnikov theory is defined in the local frame where Vc = [ p3| . In 
the context of the present theory, this is given by 



e + nUcxt + 7:[n + '2nh -n ) + — -n(v„ - v^) , 



m , 

y 



while the local energy density in the original lab frame is given by 



dab 



mh{vn 



(92) 



(93) 



Using ( pSD and (|69| ) in (|92|), one finds that the first order corrections from ri2 cancel out, leaving 

6 = ?(°) + nf/e.t + f (n^ + 2nn(") - [nf")]^) + imn(v„ - v,)^ 



(94) 



We conclude that both the local entropy density and the local energy density are determined by the diffusive local 
equilibrium distribution function /^"^ alone, and are not altered by the deviation f'-^\ 

In contrast, as we now show, the total pressure and the chemical potential must be carefully defined so that they 
satisfy the superfluid thermodynamic relations in (Bq). We first define the non-equilibrium pressure by 



P' = P+|(n2 + 2nn-n2) 



.(1) 



gric {1) 

V o ^ 12 ' 



Using ( pSD and ( |69| ) and working to first order in r-j^2 i one obtains 

P' = P 
where P is the (diffusive) local equilibrium pressure defined as 

P = P(") + ^{n^ + 2nn(") - [nW]^). 



(95) 



(96) 



(97) 



We find that the LK thermodynamic relations given in (pSl) are not satisfies if we assume that P' is the pressure (P) 
and fic is the chemical potential (/i). Extra terms appear which are associated with ri2. This means that the above 
identification of the thermodynamic variables is only valid in the lowest order hydrodynamic equations, where there 
is no dissipation. 

We recall that in deriving the Landau equations from the ZGN' equations in Ref. g], we defined the total pressure 
by (Bq) and n = /Zc, and also found extra terms in the thermodynamic relations proportional to (5/idiff (see Eq. (71) 
of Ref. 1^). Therefore the precise equivalence between the ZGN' hydrodynamics and the Landau theory shown in 
Ref. P] was restricted in the limit ujt^ — > 0, i.e., when (5/Xdiff — *■ 0. In contrast, if the pressure P is defined to be 
^ ,(0) ^ 



( P7| ) and fi = fic , ^e can show that the the superfluid thermodynamic relations in ( pq ) are satisfied. Therefore we 

conclude that the local equilibrium pressure defined in ( p7\ ) and the local equilibrium chemical potential fic given 
by (|3^) are, in fact, the correct variables to be use in the Landau-Khalatnikov equations. We will show later that 
the corrections to the total pressure and the chemical potential actually give rise to the additional damping terms 
associated with the four second viscosity coefficients Ci in (|8^). 

We now proceed to derive the LK e quations from our microscopic theory, one by one. Our continuity equatio ns fo r 
TT-c and n are given by (pa) and (67a). Adding them, we obtain the continuity equation for the total density (86a). 
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To derive the equat ion (p6b| ) for the total current j, we combine our two continuity equations and the two velocity 
equations (8b) and (|67b| ) to give 






dt dxu 



dxu 



D 



M^/ - -^Sf^^iTiD) 



(98) 



Using (R6h, we find 



5jp 



d 



dx. 



dt 



dx,, 



D^,, - -<5^.(TrD) 



+ OtifT^ o ^ 12 



(99) 



To consistently include damping due to the first order correction term in the Chapman-Enskog expansion, we use 
ri2 = T\2 as given in ([44|). We then find that ( p9| ) is identical with the LK equation (B6b) with the second viscosity 
coefficients Ci and ^2 given by 






(100) 



Using fic = fJ' 



(0) 



gn 



n(l) 



-(1) 



(1) 



H — gTT\2 (see (|35|)) and the expression for T\ y' in our equation for the condensate velocity 



given in (Eo), we find the latter can be written precisely in the LK form (S6c). Comparison between the two equations 
shows that the second viscosity coefficients ^3 and C4 are given by 



9 



gric 



Cs = — nTfiCFH, C4 = T, — Tu,aH 



3m 



(101) 



We note that our results for the second viscosities satisfy the Onsager reciprocal relation Ci = C4 (this equality follows 
quite generally, as shown by Eq. (4.28) of Ref. P3). 



Finally, we derive the equation for the local entropy. Using ( p8q ), we have 



ds C/6 un u 

dt dt ^ dt ^ ' dt^ ^ '^ 



(102) 



With the expression for the local energy density e given in (|92|), we find ( |102| ) reduces to 



ds de 



~ N ^ dn 



dt at + [^-t-M + 5(n + n)]^+5n.g^ 



dn dn m 



dt 2 



(v„ 



dt 



~,(i)^" 



dh 



dl+9--^+9-c^- 



(103) 



Here we have neglected the last term in th e first line, since it is of third order in the local velocities. Using our 
hydrodynamic equations (67), we find (|l03| ) reduces to the form ( 86d ), assuming the entropy production rate Rs is 
given by 



Rs = T-p5ri2 i 9'^ ■ K(vc - v„)] + -n^V • V, 



-2ri 



D^,, ~ :,S^,{TtD) 



^{^Tf. 



(104) 



Using ( [44|) and the expression in (|lO0| ) and (101), we see that (104) is equivalent to the Landau-Khalatnikov |12| ] 
expression given in (p7). 

We have thus shown that our equations based on a microscopic theory built on Bose condensation can be written 
in a form precisely identical to the phenomenological Landau-Khalatnikov two-fluid equations including the damping 
associated with the shear viscosity, thermal conductivity and four second viscosity coefficients. An analogous derivation 
of the Landau-Khalatnikov equations for a uniform Bose gas was first given by Kirkpatrick and Dorfman p3[ . However 
at finite temperatures, where the dominant excitations are particle-like Hartree-Fock excitations, KD did not obtain 
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the second viscosities since they neglected the source term ri2 associated with deviation from local equilibrium 
produced by the C12 collisions. We have shown that the second viscosity coefficients are directly related to the 
ri2 term first discussed in Ref. |^, which represents the collisional exchange of atoms between the condensate and 
non-condensate. 

In the above derivation of the second viscosity terms, we used ri2 — T]^2 . This restricts the validity of the results 
to the case lut^ <C f when we consider collective fluctuations with frequency to. However, our discussion can be easily 
extended to the situation when t^r^ is not small, by using (see (p^)) 

ri2(w) = , """ ( V • K(v, - v„)] + \n,V • v„| . (105) 

Using this expression, we can still write our equations in the Landau-Khalatnikov form, but now with the frequency- 
dependent second viscosity coefficients 

c,(c.) = ^: . (106) 

f - UVTfj, 

Everything else in our derivation goes through. 

The expression (ttOq) for the frequency-dependent second viscosity coefficients has in fact the expected form, as 
derived from general considerations ||25|] . The second viscosity , such as associated with compression and expansion, 
arises when a gas is coupled to an internal relaxation process (for example, the transfer of energy from the translational 
degrees of freedom of a molecule to the vibrational degrees of freedom). If the relaxation time of the internal process is 
denoted by tr, the frequency-dependent second viscosity coefficient is given by C(a;) = Co/(l — ^wtr), where Co ot tr. 
In a Bose-condensed gas at finite temperatures, the non-condensate atoms are coupled to the condensate degree of 
freedom, and we have shown that the characteristic relaxation time for the equilibration between the two components 
is given by r^. In this connection, we might recall that in the superfluid **He, the second viscosities are due to the 
fact that the total number of p hono ns an d rotons is not conserved ||T^. Above Tbec (when Uc = 0), all the second 



viscosity coefficients Ci in (fOC) and (fOf ) vanish, as expected in a normal dilute single-component gas. 

We finally note that the Landau-Khalatnikov equations could have been derived from the ZGN' hydrodynamic 
equations if we simply identified the total pressure P with that by ( ps] ) and (|9^) , and took the chemical potential to 
he n — ^c + 9T^Sri2 (see Eq. (|70|)). This leads more naturally to frequency-dependent second viscosities. On the 
other hand, the physical meaning of this choice of the local equilibrium pressure and chemical potential is not made 
clear. 

V. CALCULATION OF TRANSPORT COEFFICIENTS 

In this section, we solve the linearized equation for functions A, B and D in ( pq ) which determine the deviation 
from local equilibrium as described by V' in (Hq) and (p4). We can then calculate the transport coefficients ry and 



K as given in (62). We follow the standard procedure in the Chapman-Enskog method, as reviewed in pM- In this 
approach, one solves the linearized equation by expanding -0 in a basis set of polynomial functions. These polynomial 
functions are chosen to satisfy the constraints such as (p3f) which the solution ■0 must satisfy. In a classical gas, one 
uses Sonine polynomials [Q. One can also define analogous polynomials for a degenerate Bose gas [^. As usual, 
we calculate the transport coefficients using the lowest order polynomial approximation, which usually gives very 
accurate results for the transport coefhcients. For a more detailed mathematical discussion which is easily generalized 
to Bose condensed gases, we refer to [ p4pf| ]. 

A. The thermal conductivity 

In evaluating the thermal conductivity, it is convenient to rewrite ( |62b| ) as 

1, /" rfp ,, ^ r "^"^ 535/2(20)1 , ,_, , ,. 

Joit + ./oj 



3'^y (2^"^^")-" 



2fcBTo 253/2(2:0) 



= -\kB J J^^Mn) ■ L[nA{u)]. (107) 

Here we explicitly use the st atic equilibrium distribution function /o to evaluate the transport coefhcients. To solve 
the linear integral equation (55a) for A{u), we introduce a simple ansatz of the form |I^pl] |: 
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A{u) = A 



■mu^ 555/2(20) 



2fcBro 253/2(20) 



(108) 



This is the lowest order pol ynom ial function that satisfies the constraint given by (56a). The constant A in (lOS) is 
determined by multiplying ( ^5a| ) by u.[mv? / 2k-QT ~ 555/2(20)7253/2(20)] and integrating over p, giving 



A = 



dp 



{2T:hf 



mu^ 555/2(20) 



dp 



2fcBTo 253/2(20) 
mu^ 555/2(20 



/o(l + /o) 



{2-Khf 



2kBTo 253/2(20) 



ui 



mu'^ 555/2(20) 



Using (|l08|) and (|To|) in (|lO^, we find 



dp 



{2T:hy 



mv? 555/2(20) 



2fcBro 253/2(20) 



dp 



{2T:hf 



mv? 555/2(20) 



2fcBTo 253/2(20) 



ui 



2fcB7o 253/2(20) 
/o(l + /o) I 

TOU^ _ 555/2(20) 

2fcBro 253/2(20) 



To evaluate the p integrals in (IIC), it is convenient to introduce the dimensionless velocity variable by 



1/2 



1 



1/2 



2kBTf)J \2mkBTo^ 

With this new variable, we can rewrite the linearized collision operator as 

SmiksTofa^ 



Lm 



TT^h^ 



(^L'^M + ^''^^ricoAlL[M), 



where the dimensionless collision operators ^22 ^^^ -^12 ^^^ defined by 



^22 M = / dx2 / dx3 / dx4(5(x + X2 - X3 - X4)S{x + ^2 - ^3 - ^4) 

X/10/20(1 + /30)(1 + /40)(V'3 + -04 - V'2 - Tp) , 

^12 M = / dxi / dx2 / dx3(5(x - X2 - X3)5{xl - Pogrico -x\- 2:3) 

X [(5(x - Xi) - (5(x - X2) - (5(x - X3)](l + /io)/2o/3o(^2 + ■03 - -0), 

where /^o = {zQ^e^^ — 1)^^. Carrying out the p (or u) integrals in (|110|), one finds 



75fcB (mkBrV''^ 



64:a^T 



m \ TT 



vr 



1/2 



12^2(^0) +Ai]neo/r2(^o) 
where the functions /;f2(^o) and /f2(2o) are defined by 



'7 , , 552/2(20)' 
2 253/2(20) 



-'22 = ~ / '^-''- -"■*'' ■ ^22 [^-^ J J 



(109) 



(110) 



(111) 



(112) 



(113) 



(114) 



(115) 



(116) 



'12 



-7r^/2 / dx xa;2 • L22[xx2 



In Ref . |14[| , we derived a convenient formula for the integral /^j ? namely 

^22(^0) = \/27r^ = / dxo / dxr dy dy' F22{xo , Xr , y , y ] zq) 
Jo Jo J-l J -I 



(117) 
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xxtxl{y' + y'^-2yy% 



(118) 



where 



with 



F22{xQ,Xr,y,y']zo) 



(1 — zoe ^i)(l — zoe ^2)(1 — zqe ^3)(1 — zqe ^4)' 



(119) 



^l = 2*-^° "*" 2xoXry + x^), a;^ = -(x^ - 2xf)Xry + x^) 
2^3 = 9(^0 + 2xoXr2/' + x^), xl = -(xg - 2xoXry' + xl). 



(120) 



We note that l22{zo) is a universal function of the equihbrium fugacity zq, where zq — e Psnaoir) ^ ^q derive a similar 
expression for /fj, wc introduce the transformation 



1/ X 1/ 

X2 = 2 (^0 + Xr), X3 = -(Xq - X-r 



(121) 



We then express x^ in the polar coordinate {xr, 0, (/>) where 9 is the azimuthal angle with respect to the vector Xq, ie, 
Xr • Xq = XrXo COS 9 . With these new variables, one obtains the following expression for /fj: 



7^2(^0) = 87r^/2 / dxr 
Jo 



dyFi2{xr, y; zo)xl{xl + Pogucaf/'^ 



xl{xl + 3/3osnco)(l - y^) + t(/3o5"co)^ 



where y = cos 9 and 



with 



Fl2{Xr,y\Zo) 



ZqC 



(1-zoe ^i)(l-zoe ^2)(1 - zoe =^3)' 



(122) 



(123) 



xl = 2{xl + Pogrico), 



xj=xl + Xry\/x^ + (iogrico + -zPogncO 
xl^xl- Xry\/xl + PogricQ + -Pognco- 



(124) 



The formula in (115) gives the thermal conductivity k as a universal function of gnco(r) or equivalently in ter ms o f 
the local fugacity zo{r) = e~^°^^''°^'^\ If we ignore the contribution from C12 collisions, i.e., set /fj to zero, (115) 
reduces to the expression for k derived in our earl ier work |l4,nffl. 

One can also write the expression for k in (|115[) in the following useful form: 



5V2 (^ nnklTo \ j 7gr/2izo) 
m J \ 2^5/2(20) 



5^5/2(^:0) 
2.93/2(2:0) 



(125) 



where r^ is the "thermal relaxation time" associated with the thermal conductivity, as defined in Appendix ^. In turn, 
one can also write the reciprocal of this relaxation time Tk as the sum of contributions from C12 and C22 collisions, 



1 



1 



1 



7k, 12 Tk.22 



(126) 



where these relaxation times are given explicitly in Appendix ^. The physical meaning of this r^ relaxation time is 
discussed in Appendix H, using a simple relaxation time approximation for the collision integrals in (11]) . 
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B. The shear viscosity 



In evaluating the shear viscosity 77, it is convenient to rewrite (62a) as 



m 



10 / (2^^ yuf^u^ - -S^^u^j B{u) iu^u^, - -(5^^u2 ) /o(l + /o) 



fee To 



dp 



10 J (27r?i)3 



B{u) ( Uf,u^ - -Sf^tyu"^ ) L 



B{u) I u^u^ ~ -;jS^„u 



(127) 



To solve (127), the simplest consistent ansatz |l4|2l| is to use B{u) = B. The constant B can be determined by 
multiplying (55b) by [u^j^Uu — <5^i/U^/3) and integrating over p, to give 



B 



m 



dp 



^bTo y (27ra)3 V " "^ 3 



Uf^Uy ~ ^S^ivV^ ] /o(l + /o) 



(2^ (-.-. - \^,^u' ) L 



u^u^ - -S^„u'^ 



(128) 



Using this in (127), we obtain 



V^ 



dp 



IO/cbTo W {2Trh) 



dp 



3 \ "■t' 



Uf^u^ - -Sfj.i,u-' 1 /o(l + /o) 



(27r?i)3 V '' 3 ^" 



(5„,,M L 



Ufj^Ui, - -S^jyu'^ 



With the dimcnsionless variable defined in (11_1), this expression for the shear viscosity 77 can be rewritten as 



'^ 32V2a2^ "^ °^ /2''2(zo)+A3n,o/I'2(^o) 



(129) 



(130) 



where 



-^22 — / '^■^ 1 '^i^-^i' o M^'^ ) -^22 



Xl_iXl/ ^0 jj^yX 



3 /"OO /"OC /"l pi 

-j= j dxo dxr dy dy' Fi2{xQ,Xr,y,y'] zq) 
xxlxl{l + y^+y'^-iy^y'^), 



(131) 



and 



-^12 = -7^^^^ / cix f X^,Xi, ~ -5^,^x^ j L'^2 

nOO nl 

= Stt"^/^ / da;^ / dyFi2{xr,y;zo) 
Jo J-l 






xoi^a/x^ +/3ognco 



x^(x^ + /?off«co)(l -y ) + g(/3o5«co) 



(132) 



These expression s inv olve the same functions F22 and F12 defined earlier in (11£) and ( Il23|). 

Analogous to (125) for the thermal conductivity, can also write the expression for 77 in ( 130 ) in the following form: 

"55/2(20)" 



77 = T^hokBTo 



53/2(20) 



(133) 



where the viscous relaxation time r^ is defined in Appendix H. Again one can write the reciprocal of the relaxation 



time T^i as 



'''■ij,12 'Tri.22 



(134) 



where these C12 and C22 relaxation times are given in Appendix H. As with r^, the relaxation time r^ can be 
understood in terms of a simple relaxation time approximation for the collision term in the kinetic equation. 
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C. The second viscosity coefficients 



To find the expression for t as defined in (b3), we use the simple ansatz for the form of the solution for D{u) of 

(Hi), 



Diu) = £» (72 



3 fee To 



CTl 



(135) 



where cti and a2 are defined in (5l[ ). As usual, we leave the dependence on (r,t) implicit. One easily verifies that 
(135) satisfies the constraint (B6b). The constant D can be determined by integrating (55c) over p: 



D 



dp 



{2'Kny 



:Li2[cT2 



mu 



■<Ji] 



1 



dp 



n(") J (27r?i)3 
dp 



(72 



3 fee To 



/o(l + /o) 



'Pg9t 



rii2[l]. 



{2nnf 
Using Li2[^] = /3o5"'coT/i2[l] and / |2^E^ii2[l] = -nco/n2, where ri2 is defined in (|2|), ( p^ ) gives finally 

(ti2 + PogricoT) 



(27Thy^ 
D 



ncQ[(J2 + |o-i/3(/Xco - f/o)] 



(136) 



(137) 



Using (137) and (|l35|) in the expression for r given by (63), we can solve to give an explicit expression for r, namely 



£ _ J_ r nco 

T Tl2 \ no 



0-2 + -cri/3o(McQ - ^o) 



Pogric 



<7HTi2 



(138) 



We thus see that r is precisely the relaxation time t^ first introduced in the ZGN' two-fluid hydrodynamics. We can 
now express the four second viscosity coefficients in (100) and (101) in terms of the T12 collision time defined in (pj 



Ci 



fceT 2 /■ n-ckBT ^ /- 

-5 — o-^n2, C2 = — 7. — cr^n2, Cs 
3m 9 



mPric 



<^h'''i2, C4 = Ci- 



(139) 



D. Numerical results for a uniform Bose gas 



For illustration, we calculate the transport coefficients in (lllq) and (130) for a uniform Bose gas (C/oxt =0). As in 
Ref. Q, we choose gn/fceTBEC = 0.2. In Fig.l and Fig. 2, we plot the temperature dependence of the dimensionless 
transport coefficients R and fj, defined by 



n/nv^iTokB, rj = T]/nVcimTo. 



(140) 



Here Vd = (5fcBTBEc/3TO)^^^ is the hydrodynamic sound velocity of a classical gas at T = Tbec and Tq"^ = 
-\/2(87ra^)n(8fcBrBEc/7''m)^/^ is the classical mean collision time evaluated at T = Tbec- To see the separate ef- 
fects of the C22 collisions and C12 collisions, we also plot the results obtained by taking either 1^2^ = or /^j'' = 0- I^^ 
Ref. pOl, we neglected C12. We see that the both k and 77 are reduced when we include the C12 collisions. At lower 
temperatures T ^ 0.5Tbec, both k and 77 are dominated by the C12 collision integral. 

In Fig. 3, we plot the four second viscosity coefficients given in (13£) for a uniform Bose gas. We recall that in 
Fig.l of Ref. |g[, we gave the temperature dependence of r^ and T12 for gn/kBTBEC = 0.1. In Fig. 3, we use the 
dimensionless second viscosity coefficients, defined by 



Ci = Ci/WciT-Q, C2 = C,2m/nvlyTo, Cs = C2n-/™WciTo- 



(141) 



The transport coefficients in a trapped Bose gas behave quite differently from those of a uniform Bose gas. In 
particular, since Uc^ h holds in the central region of the trap at all temperatures, the contribution of C12 collisions 
dominates over the contribution of the C22 collisions at all temperatures below Tbec- We will discuss the implications 
of this at the end of Section VI . 



20 



VI. CONCLUSIONS 

In this paper, we have derived two-fluid hydrodynamic equations starting from the quantum kinetic equation and 
the generahzed GP equation derived in |g,g,|lj. However, to complement and extend our earher work |||,g|, we started 
from the complete local equilibrium single-particle distribution /*•"■' as given by ( |37| ) and (Bq). Using the Chapman- 
Enskog approach, we then included the effects of a small deviation from this local equilibrium form. This deviation 
from local equilibrium within the thermal cloud brings in the usual kind of hydrodynamic damping due to the thermal 
conductivity and shear viscosity of the thermal cloud. A summary of our major results is given in the final paragraph 
of Section |. 

In addition, we showed that the there is additional dissipation associated with the coUisional exchange of atoms 
between the condensate and non-condensate components. When we write our hydrodynamic equations in the Landau- 



Khalatnikov 1 12 2^ form given by (|8^) , this damping is described in terms of the usual four second viscosity coefficients 



for a Bose superfluid. The appearance of the second viscosity coefficients in the equations for the total current j in 



( p6q ) and for the superfluid velocity Vc in ( B6c ) is due to the deviation of the total pressure and the chemical potential 
from their local equilibrium values. We might also recall that Khalatnikov ||l^ discusses a specific model for the second 
viscosity coefficients in superfluid ^He by introducing "local chemical potentials" for the phonons (/iph) and rotons 
(/ir). These describe a situation where such excitations (describing the normal fluid) are out of local equilibrium with 
the superfluid component. Clearly this discussion has connections with our calculations based on the condensate and 
non-condensate not being in diffusive equilibrium. 

The frequency-dependence of the second viscosity coefficients is a result of the fact that our two- fluid hydrodynamics 
deals with the dynamics of the condensate and non-condensate components as separate degrees of freedom. This feature 
is made more explicit in our recent papers ||3M|lO|. In particular, it gives rise to a new relaxational zero frequency 
mode. As mentioned at the end of Section y, and more explicitly in Ref. [|0|, this mode may be viewed as the 
(renormalized) version of the zero- frequency thermal diffusion mode |20[| above Tbec- The presence of this new mode 
below Tbec is somewhat hidden in the formulation in terms of the LK two-fluid equations given in ( pq ) . 

In Section 0, we derived explicit formulas for all the transport coefficients within our model. In Ref. [Q, we 
only took into account the deviations from local equilibrium due to the C22 collision integral in calculating the shear 
viscosity and the thermal conductivity coefficients. In the present paper, we have also included the contribution to 



these quantities from the C12 collision integral. From (125) and (133), we see that both k and ry are given in a form 



proportional to characteristic relaxation times Tk = (t^Tj^j + t^22)^"^ ^^"^ '''n — (''')7i2 + ''^sa)" ' respectively, which 
are defined and motivated in Appendix g. 

In a rough estimate, w e fin d t„^22 ,''>), 22 ~ ^ci ~ ^/n and Tk,i2,t,,^i2 ~ {n/nc)Tci ~ l/«c, where Tc\ is the classical 



collision time defined in (|B2[) . We therefore observe that the effect of C12 collisions reduces the magnitude of both 
K and ?7 by a factor ~ 1/(1 + ric/n), a result also noted in Ref. [ p3[ for a uniform gas. The contribution of the C12 
collisions is always important in a trapped gas, since the condensate density at the central region of a trap is much 
larger than the density of the thermal cloud even at temperatures very close to Tbec- In a trapped Bose gas, we 
find the effect of the C12 collisions is enhanced by a large factor ric/n ^ 1. This means that k and ry are always 
dominated by the contribution of the C12 collisions. Since the effective relaxation times Ti^^i2 and t^,i2 are smaller 
than the classical collision time t^i by a factor n/ric <^ 1, this implies that in a trapped Bose gas, the hydrodynamic 
region may be much easier to reach at finite temperatures than expected from naive considerations based on using 
the classical collision time (i.e., toTd <ti 1). That is to say, one might easily have u;Tk,i2 ^ 1 and wt^,i2 "C 1, even 
though LOTci 3> 1. This has very important implications in deciding if one is in the coUisionless or the hydrodynamic 
region. 

One problem not dealt with in this paper is the fact that in a trapped Bose gas, the decreasing density in the 
tail of the thermal cloud means that the hydrodynamic description breakdowns eventually. This problem enters the 
evaluation of the expressions for the rj and k transport coefficients given in Section [V|. In recent papers dealing with 
the case above Tbec p7| , [L4[ , this problem was handled in an ad-hoc manner by introducing a spatial cutoff in the 
thermal cloud. In a future paper, we given an alternative approach based on starting with an improved solution of 
the kinetic equation, which naturally includes the cross-over to the non-interacting gas in the thermal gas tail. 
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APPENDIX A: 



We briefly sketch the derivation of the first-order kinetic equation given in (M. Using (pTl) in (M3), one has 



h — • Vr 

dt m 



Vrt/(r,t)-Vp 



/^°Hr,p,t) 



z\dt^ m' )^ ' 2fcBT2 \dt ' m 



? + i^.V|T 



mw. 
kWT 



— + -^ • V I v„ 
ot m 



Vt/(r,t) 
knT 



r>ii+r>). 



(o)n 



(Al) 



The notation d'^ /dt is explained below (|4S|). Using the expressions for the density ri'^) given by ( [40| ) and the pressure 
P(0) in (H), one finds 

dt ^ 2T "ar^~ aT' 



9t 



2r 9i 



9t ' 



(A2) 



where 7 is the v aria ble introduced after (bW and z — z^^"^ as defined below (Eoh. One may combine these equations 
with (43a) and (43c) to show that the equations in (A2) reduce to 



'df 

d°z 



;T(v-v„)-v„-vr- 



T7 I ^12 

-v„ • Vz + a2Z-—yr 



2T 

3n(") 



CTiri 



(A3) 



The analogous equation for d'^Vn/dt is given directly by (43b). Using these results in (Al), one finds that it reduces 
to 



cP_ p 

dt rr, 

T 

0'2 



Vr(7-Vp / 



p(o) 



VT 



5P(0) 



2kBT 2n(o)fcBT / /cbT 



SfceT 



Cl 



fcRT 



U • (Vc - V„) 






u- (u- V)v„- —V • v„ 

/(0)(i + y(0))^ 



(A4) 



where we recall u = p/m — v„. In calculating the dissipative terms, we only consider terms to first order in the 
velocity fields v„ and Vc. Since ri2 is proportional to v„ and Vc (see (p4|)), we can neglect the last term (proportional 
to Vc — v„) in (A4). This linearized version of ( A4) can be rewritten in the form shown on the left hand side of (pQ). 
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The relaxation time r^ in (125) is defined by 



dp 



(27r?i)3 
dp 

{2TTn)^ 

157r9/2a' 
8TO(fcBro)2a2 



mu2 555/2(^0) 



2fcBTo 253/2(^0) 
mu2 535/2(2:0) 



«'/o(l + /o) 



uL 



'( rmP 535/2(2:0) 



2fcBTo 2^3/2(20) 

1^2(^0) + ricoAlIuizo) 



\2fcBr0 2^3/2(20) 
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15\/2^^/2 



-Tcl 



Here 



297/2(^0)33/2(^0) - 2.95/2(^0) 



\/2(87ra2)no(r)(8fcBro/7rm)i/^ 



(Bl) 



(B2) 



is the collision time of a classical gas with density fiQ and quantum cross section a = Stto^. In turn, the relaxation 
times in (|126|) are defined by 



Tk.22 



dp 



and 



TkA2 



dp 

{2TThf 

8m(fcB7o)2a2 

15V2^V2 
: '^ci 



dp 



mu^ 5^5/2(20) 



2fcBTo 2^3/2(20) 

mu^ _ 5g5/2(^o) 
2fcBTo 2^3/2(^0) 



«'/o(l + /o) 



u • L-} 



mv^ 555/2(2:0) 



2fcBTo 2^3/2(^0) 



57/2(^0) - f35/2(^o)/ff3/2(2;o) 



^2^2(^0) 
57/2(^0)33/2(20) - 1.95/2(^0) 



(27r;i)3 
dp 



^22(^0) 



mu^ 535/2(2:0) 



2fcBT'o 2^3/2(20) 

mv^ _ 5g5/2(2o) 

{2^nf L2fcBTo ~ 2^3/2(20) 



^i'/o(l + /o) 



U-ii2 



mu^ 555/2(20) 



2fcBTo 2^3/2(20) 



" 8TO(fcBTo)2a2 
In an analogous way, the relaxation time r,, in (|133|) is defined by 



2^7/2(20) ~ ig5/2(2o)/g3/2(2o) 

ncoA3/f2(2o) 
57/2(20)93/2(20) - §55/2(^0) 



ncoAi^/f2(2o) 



3 ( U^^Ui, - -S^^u^ ) /o(f + /o) 



{2t:K) 
dp 



2^^)3 1 "'^"'^ ~ g^A"'" ) ^ 



I 



Uf^Uu - -S^uW 



5^9/2^3 



16m(/cB2o)2a2 

5V27r7/2 



55/2(20) 



-Tcl 



^22(20) +ncoAg/l''2 (20) 
55/2(20)53/2(20) 



/2''2(2o)+ncoAj]/i''2(2o) 



The relaxation times in (134) are defined by 

dp 

T^,22 



(27rS.)3 * "^"^ ^ ^'^^'''"^ ) •^°^^ + -^^^ 



(B3) 



(B4) 



(B5) 
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Ufj.u„ ~ -Sfj_^i 



(B6) 



and 



T,,,12 



■^ [Uf,u^- -S^^u'^ ] foil + fa) 



{2TTh) 

dp 



{2TTnf 

16m(/cB2o)^a^ 



- Ufj^Uiy ~ i^S^i^u 1 L 



-Sui^U 



95/2{zo) 
UcoAll^'^izo) 



5V27r7/2 



-Tcl 



g5/2{zo)g3/2{zo) 

ncoAg/{'2(zo) 
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(B7) 



Tcl (no (r)) and 



We note that in a non-degenerate gas, these expressions simphfy and we find r^ — t^^22 
Trj — TV;, 22 — tTci(no(r)). The latter expression for r^ agree with the result for the shear viscosity given in Eq. (11) 
by Kavoulakis et al. p8[ |. 

The physical meaning of these new relaxation times becomes clear when we compare our Chapman-Enskog analysis 
in the text with a simple relaxation time approximation |20[ . In the relaxation time approximation, one simply 
replaces the collision term in (^ with — [/ — /^"•'I/tioi, where TioI is a phenomenological relaxation time characterizing 
ho w fas t the syste m relaxes to local equilibrium. With this approximation, the solutions of the linearized equations 
in ( 55a ) and ( p5b| ) for the functions A{u) and B{u) are found to be simply given by 



A{u) 



-Trcl 



mv? 5.95/2(2) 



2fcBT 253/2(2) 



B{u) 



-Trcl 



fceT 



(B8) 



In contrast, our Chapman-Enskog solution for A{u ) is given by (108) with the coefficient A given by (109), while one 
has B (u) — B with the constant B given by (|12 
( |108| ) and (12?) can be written as 



In terms the relaxation times t^ and r,, defined above, we find 



A{u) 



mu^ 5g5/2{z) 



2k^T 2^3/2(2) 



B{u) 



'k^T- 



(B9) 



Comparing (B9) with (B8), we see that both r„ and t^ can be identified with the relaxation time Tici- That is, in the 



simple relaxation time approximation, k and 77 are still given by the formulas (125) and ( |133| ), but with r^ — r,, = Trci 



We also note that Eq. (7) of Rcf. |^ gives a general expression for various coUisional relaxation times above TbeCj 
which is given by a formula analogous to (Bl) and (B5). 

It is these effective relaxation times that determine whether one is in the hydrodynamic two-fluid domain, lot,^ < 1, 
^Tii ^ 1 (see discussion in Section VI). They should be used in place of the classical colhsion time Td in (B2). 
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Hohenberg and Martin p3| worked out the dispersion relation of the hydrodynamic modes in a a uniform Bose 
superfiuid using the Landau-Khalatnikov two-fluid equations. The frequencies of the first and second sound modes 
are given by 



J = Uik — iDjk Lo, 



(CI) 
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where the sound velocities Ui are determined by the coupled equations 

Tpss^ dP 



2 I 2 



2 2 
U1U2 



PnCv ds 

Tp,s^ dP 

PnCv dp 



(C2) 
(C3) 



Here s = s/mn is the entropy per unit mass and Cy is the specific heat per unit mass. The damping coefficients Di 
are determined by the coupled equations 



Di + D2 = h (Ci + C4) + P + -^: 

opn Pn Pn Pn PCy 



ulDi + ulD2 



dP C2 + 477/3 ps /Ts2 2Ts dP 
dp p Pn \ Cy pcy dT 



dp 

dp 



+ C; 



P^dP 

Pn dp 



Pn \ Cyp dT 



dP 
dp 



(C4) 
(C5) 

(C6) 



We note that the above general expressions are valid for both liquid ''He and Bose gases. For the liquid ■*He, these 



formulas can be simplified by using Cy p^ Cp. One has 

2 dP 

2 ^, P^ dT 



D2 



Pn 9(l/s)' 

C2 + 47?/3 

p 

1 r K 



Ps 



p [Tids/dT) Pn 



Cap' - p(Ci + C4) + C2 + x?? 



(C7) 

(C8) 

(C9) 

(CIO) 



However, these are not valid for a dilute Bose gas. Calculating the various thermodynamic derivatives with the 
Hartree-Fock single-particle spectrum used in this paper ||26| , the sound velocities and damping coefficients for a 
uniform Bose gas are given by 



ul 



u\ 



Di 



5Po , 2(7^0 [^ 2n2 



gn-co 



m 

{i-<yH) 



^"^ 



m 
:^ + ^-^(Ci+C4) + ^ 

oPn Pn Pn Pn 



4 kTo 



9 PnUJ 



2ric0 
3no 



P>2 ~ PsCa + 7:-^i(^i'^H)^KTo 



Here an is defined in 



Qplui 

while Pn — mfiQ and ps — mrico. 



(Cll) 
(C12) 

(C13) 
(C14) 
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FIGURE CAPTIONS 



FIG.l: Plot of the thermal conductivity k in a uniform gas for gn = 0.2fcBTBEC as a function of temperature. We also 
plot the results by taking /fj = (dashed line) and I22 = (dot-dashed line). 

FIG. 2: Plot of the shear viscosity coefficient in a uniform gas for gn — Q.2kBTBEC as a function of temperature. We 
also plot the results by taking 1^2 = (dashed line) and I22 — (dot-dashed line). 

FIG. 3: Plot of the second viscosity coefficients Q in a uniform gas for gn = 0.2fcBTBEC as a function of temperature. 
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